Scaling analysis of Schottky barriers at metal-embedded semiconducting carbon 

nanotube interfaces 
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We present an atomistic self-consistent tight-binding study of the electronic and transport prop- 
erties of metal-semiconducting carbon nanotube interfaces as a function of the nanotube channel 
length when the end of the nanotube wire is buried inside the electrodes. We show that the lineup of 
the nanotube band structure relative to the metal Fermi-level depends strongly on the metal work 
function but weakly on the details of the interface. We analyze the length-dependent transport 
characteristics, which predicts a transition from tunneling to thermally-activated transport with 
increasing nanotube channel length. 

PACS numbers: 73.63.-b,73.40.-c,85.65.-|-h 
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The nature of Schottky barrier formation and its ef- 
fect on charge transport through metal- semiconductor 
interfaces have been actively investigated for decades 
due to their importance in microelectronics technology, 
but are not fully resolved despite the tremendous ef- 
forts from both the experimental and theoretical sides i 
The rapid development of single-wall carbon nanotubes 
(SWNT) as a promising device technology functioning 
at the nano/molecular-scale^-'^-'' presents a new chal- 
lengei^'S'i SWNTs are nanometer-diameter all-carbon 
cylinders. Unlike the planar metal-semiconductor inter- 
face, both the contact area and the active device region in 
a metal-nanotube (NT) interface can have atomic-scale 
dimensions. In addition, due to the weaker effects of 
electron-impurity* and electron-phonon scattering^' in 
quasi-one-dimensional systems, transport through nan- 
otube junctions can be either coherent or phonon-limited 
depending on the nanotube type, channel length, temper- 
ature and the bias voltage. These have made it difficult to 
assess the Schottky barrier effect on the measured trans- 
port characteristics. 

In this paper, we present an analysis of Schottky bar- 
rier formation at a model metal-SWNT interface as a 
function of the SWNT channel length. The model sys- 
tem is illustrated schematically in Fig.Q where the ends 
of an infinitely long SWNT wire are buried inside the 
semi-inifinite metallic electrodes and the channel length 
is determined by the distance between the source and 
drain electrodes JS- The embedded contact scheme is fa- 
vorable for the formation of low-resistance contacti^i^tiS 
We choose (10, 0) SWNT as the protype semiconduct- 
ing SWNT, whose work function is taken as that of 
the graphite (4.5 eV)i^ The SWNT channel length in- 
vesitigated ranges from L = 2.0, 4.1, 8.4, 12.6, 16. 9(nm) 
to 21.2 (nm), corresponding to number of unitcells of 
5,10,20,30,40 and 50 respectively. We calculate the 
transport characteristics within the coherent transport 
regime, as appropriate for such short nanotubes We 
consider gold (Au) and titanium (Ti) electrodes as ex- 
amples of high- and low- work function metals (5.1 and 
4.33 eV respectively for polycrystalline materials). In 
view of recent report on ballistic transport in palladium 



(Pd) contacted SWNT devices^ we also consider Pd 
electrodes with a similar workfunction to gold (5.12 eV 
for polycrystalline materials) liSi 

We analyze the Schottky barrier effect at the metal- 
SWNT interface by examining the electrostatics, the 
band lineup and the conductance of the metal-SWNT 
wire-metal junction as a function of the SWNT chan- 
nel length. This is investigated using a Green's func- 
tion based self-consistent tight-binding (SCTB) theory, 
which takes fully into account the three-dimensional elec- 
trostatics and the atomic-scale electronic structure of 
the SWNT junctions. The SCTB model is essentially 
the semiempirical implementation of the self-consistent 
Matrix Green's function (SCMGF) method for ah ini- 
tio modeling of molecular-scale devices-'^^ and has been 
used previously to investigate Schottky barrier forma- 
tion when a finite-size SWNT molecule is contacted to 
the electrodes through the dangling bonds at the endii^ 

The SCTB model starts with the semi-emprical Hamil- 
tonian Hq of the bare (10,0) SWNT wire using the 
Extended Huckel Theory (EHT) with non-orthogonal 
{sp) basis sets 0m(r),^^ which gives a band gap of 
K, Q.9{eV)}^ We describe the interaction between the 
SWNT channel and the rest of the junction using matrix 
self-energy operators and calculate the density matrix 
Pij and therefore the electron density of the equilibrium 
SWNT junction from 

G« = {{E + iQ+)S-H-^l{E)~^l:Nt{E)-^r{E)-^i 
dE 
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Here S is overlap matrix and f{E — Ep) is the Fermi 
distribution in the electrodes. The matrix self-energy 
operators include contributions from coupling both with 
the left (right) electrodes 'Sl{r) ^nd with the part of the 
SWNT wire embedded within the left (right) electrode 
'^L{R);NT: which are calculated from the surface Green's 
function of the semi-infinite metal surfaces and semi- 
inifinite SWNT wires respectivelyi^^i^*':^* The Hamilto- 
nian of the SWNT channel is now H = Hq + 6V[Sp] 
where 6p is the density of transferred charge and SV is 
the induced change in the electrostatic potential. 
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FIG. 1: (Color online) Schematic illustration of the metal- 
SWNT wire-metal junction. The ends of the long SWNT wire 
are surrounded entirely by the semi-infinite electrodes, with 
only a finite segment being sandwiched between the electrodes 
(defined as the channel). Also shown is the coordinate system 
of the nanotube junction. 



The self-consistent calculation proceeds by approxi- 
mating the charge distribution as superposition of atom- 
centered charge distributions^^ Sp{r) = 5NiPi(r — ri), 
where 5Ni — {pS)ii—Nf and is the number of valence 
electrons on atomic-site i of the bare SWNT wire. Piir) is 
a normalized Slater-type function, iiiiS whose exponent 
is chosen such that / drdr' pi{r)pi{r')/\r — r'\ = li— 
where Ii{Ai) are the atomic electron affinity (ioniza- 
tion potential). This gives 5V{f) — '^^5NiVi{r — Ti), 
where Vi — J df" pi{r' — ri)/\r — r'\ can be evaluated an- 
alytically. ^'''^^ We take into account the image-potential 
effect by including within 5V contributions from both 
atom-centered charges and their image charges (centered 
around the image positions), rather than imposing an 
image-type potential correction. The self-consistent cy- 
cle is completed by calculating the matrix elements of the 
potential 5Vmn — j d'r(t)*m{'i^5V [r)(j)n{'r) using two types 
of scheme: (1) If m,n belong to the same atomic site i, 
we calculate it by direct numerical integration; (2) if m, n 
belong to different atomic sites, we use the approxima- 
tion 5Vmn = Smni^V-mm. + (SKm)/2. 

The calculated charge transfer and electrostatic po- 
tential change along the cylindrical surface of the SWNT 
for the Au/Ti/Pd-SWNT-Au/Ti/Pd junctions are plot- 
ted in Fig. (|2Jl. The electrostatic potential change is the 
difference between the electrostatic potentials within the 
SWNT junction and the bare SWNT, obtained as the su- 
perposition of contributions from the transfered charges 
(plus their image charges) throughout the junction. Here 
it is important to separate the electronic processes at the 
interface and inside the channel. The coupling with elec- 
trodes in the embedded part of the SWNT wire induces 
only a localized perturbation to the SWNT channel sand- 
wiched between the electrodes, so the electronic states in 



-0.1 



t5-0.1 



3.5 



Ti-SWNT-Ti 
X 10"^ 



10 11 



— — — j| — ^ 



10 



15 



20 




o 
a. 

E 
o 

rS 

£0.02 



(c) 



10 

Pd-SWNT-Pd 



15 



20 



-5.5 




-6 







I °HH — \ '° — ^1 — i 



■830.02^ 



10 



15 



20 




5 10 15 

Position along NT axis (nm) 



FIG. 2: Electrostatics of the three metal-SWNT-metal junc- 
tions as a function of SWNT channel length for six different 
lengths. For each junction, the inset shows the magnified 
view of the transfered charge in the middle of the channel 
for the longest SWNT studied. The electrostatic potential 
change shown here is that taken at the cylindrical surface of 
the SWNT. The horizontal lines in the potential plot denote 
the work function differences between the electrodes and the 
SWNT wire. 



the middle of the channel are essentially identical to those 
of the bulk (infinitely long) SWNT except for the shortest 
(2.1 nm) channel length studied here, leading to similar 
charge transfer both at the interface and in the middle 
of the channel. For such SWNT channels (longer than 
4.2(nm)), the potential change in the middle and conse- 
quently the band lineup scheme become independent of 
the channel lengths. Note that the transfered charge in 
the middle of the channel shows oscillatory behavior due 
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FIG. 3: (Color online) Fig. (a) shows cross sectional view 
of electrostatic potential change at the Au-SWNT wire-Au 
junction for SWNT channel length of 8.4(nm). The SWNT 
diameter is 0.8(nm). (b) shows the corresponding local den- 
sity of states (LDOS) as a function of position along the NT 
axis for SWNT channel length of 8.4(nm). The plotted LDOS 
is obtained by summing over the 10 atoms of each carbon ring 
of the (10, 0) SWNT. Note that each cut along the energy axis 
at a given axial position gives the LDOS of the corresponding 
carbon ring and each cut along the position axis at a given 
energy gives the corresponing band shift. 



to the two-sublattice structure of zigzag tubesii^ 

Compared with the finite SWNT molecule contacted to 
the electrode surfaces through the ring of dangling-bond 
end atomSfi^ the magnitude of charge transfer at the 
interface is smaller for both the Au-SWNT-Au and Ti- 
SWNT-Ti junctions since no dangling bonds are involved 
in the embedded contact scheme. But the magnitudes 
of both the transfered charge and electrostatic potential 
change in the middle of the SWNT are close to those ob- 
tained in the end contact scheme once the finite SWNT 
molecule has reached the bulk limit (20 nm and longer) 
In addition, although the transfered charge and electro- 
static potential change are larger at the interface for the 
Pd-SWNT-Pd junction due to the more directional d or- 
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FIG. 4: Room temperature conductance of the metal-SWNT 
wire-metal junction as a function of SWNT channel length. 



bitals of Pd surface atoms, their values in the middle 
of the channel are similar to those of the Au-SWNT-Au 
junction due to the similar workfunctions of Au and Pd. 
For the high (low) workfunction metals as Au/Pd (Ti), 
there is a small decrease (increase) of electron density in 
the SWNT channel, which can be identified as contact- 
induced hole (electron) doping. Therefore, the lineup of 
the SWNT band in the middle of the channel depends 
strongly on the metal work functions, but only weakly 
on the nature of the contact and the interaction across 
the interface. 

Leonard and Tersoff& first pointed out correctly that 
the strength of the interface coupling doesn't affect the 
lineup of the SWNT band in regions moving away from 
the interface, based on the electrostatics of an ideal cylin- 
der and the confined geometry of the metal-nanotube in- 
terface. However, the electrostatics of any nanostruture 
is three-dimensional. For SWNTs the electrostatic po- 
tential change induced by the charge transfer varies in 
directions both along and perperdicular to the SWNT 
axis {Z axis). This is clearly seen in Fig. |2Ia), where we 
show a cross sectional view of the electrostatic potential 
change in the XZ plane. For the (10, 0) SWNT with 
a diameter of « 0.8(nm), the change in the electrostatic 
potential inside such narrow cylinder is small, but decays 
to about 1 /4 of its value at the cylindrical center 1 (nm) 
away from the SWNT surface. 

The confined cylindrical geometry and three- 
dimensional electrostatics of the metal-SWNT interface 
lead to a completely different physical picture of the 
band shift from that of the planar metal-semiconductor 
interface. In particular, the shift of the SWNT band edge 
in the direction along the SWNT axis doesn't follow the 
change in the electrostatic potential. This is illustrated in 
the local density of states (LDOS) of the Au-SWNT-Au 
junction as a function of position along the SWNT axis 
in Fig. inib). Note that despite a w O.A{V) change of 
electrostatic potential within 3(nm) of the interface. 
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both the conduction band and valence band edges are 
nearly contant along the NT axisi^S From the LDOS in 
the middle of the channel, we can determine that for 
both the Au-SWNT-Au and Pd-SWNT-Pd junction the 
Fermi-level is located slightly below (by « 0.05(ey)) the 
midgap, while for Ti-SWNT-Ti junction it is located 
above (by « 0.15{eV)) the midgap. For the case of Au 
and Ti electrodes, these values are essentially identical 
to those obtained in the end contact scheme>i^ 

The physical principles of Schottky barrier formation 
at the metal-SWNT interface can be understood as fol- 
lows: Since the electrochemical potential (Fermi-level) 
and therefore the electron occupation in the junction 
are determined by the electrodes in contact with the 
SWNT,2i the band lineup inside the SWNT channel is 
determined by the self-consistent charge transfer effect 
through the entire metal-SWNT-metal junction. There- 
fore, the metal Fermi-level position should be close to 
the middle of the gap since otherwise extensive charge 
transfer will occur inside the SWNT channel. Since 
the screening of the work function difference inside the 
SWNT junction is weak, the metal Fermi-level should 
be below (above) the middle of the gap for a high (low) 
workfunction metal so that the net decrease (increase) of 
electrons inside the channel shifts the SWNT band edge 
down (up) relative to the metal Fermi-level. Exactly how 
this is achieved from the interface to the middle of the 
channel will depend on the details of the contact without 
affecting the lineup scheme in the bulk region. 

Given the potential shift across the metal- 
SWNT interface, we can evaluate the length 
and temperature dependence of the SWNT junc- 
tion conductance using the Landauer formula 

G - ^ J dETiE)[-^iE - Ef)] = Gtu + Gth 
and T{E) = Tr[TLiE)G^{E)TB.iE)G^iE)]}^-^^ Here 
we have separated the conductance into tunneling 
contribution Gtu — ^j;-T{Ep) and thermal-activation 
contribution Gth = G ~ Gtu- The result at room 
temperature is shown in Fig. 0] The tunneling con- 
ductance for all three junctions decreases exponentially 
with the SWNT channel length for channel lengths 



longer than 2.1(nm). A separation of contact and bulk 
effect on the tunneling resistance can thus be achieved 
using R = Rqc'^^, where Rq is the contact resistance 
and d is the inverse tunneling decay length. We find 
Ro = 8.2, 0.68, 12.3(fcf7) and d = 1.68, 1.37, 1.68(l/nm) 
for the Au-SWNT-Au, Ti-SWNT-Ti and Pd-SWNT-Pd 
junctions respectively. Note that the contact resistance is 
different for the Au and Pd electrodes due to the different 
interface coupling, but the inverse decay length (a bulk- 
related parameter) is the same. The room-temperature 
conductance saturates with increasing SWNT length, 
because tunneling is exponentially suppressed and 
transport becomes dominated by thermal-activation 
over the top of the potential barrier, whose height 
is approximately independent of the SWNT channel 
length. For Ti-SWNT-Ti junction, this leads to a tran- 
sition from tunneling to thermally-activated transport 
at roughly 4(nm). For Au/Pd-SWNT-Au/Pd junctions, 
this transition occurs at channel length of roughly 
9{nm). Consequently we can say that transport through 
the SWNT junction is bulk-limited at low temperature, 
but is contact-limited at room temperature. 

In conclusion, we have presented an atomistic real- 
space analysis of Schottky barrier formation at metal- 
SWNT interfaces with embedded contact. Further anal- 
ysis is needed that treat both the gate and source/drain 
field self-consistently within the SWNT junction, to 
achieve a thorough understanding of nanotube-based de- 
vices. 
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